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Motivated by recent advances in fabricating graphene nanostructures, we find that an electron 
can be trapped in ^-shaped graphene nanoconstriction with zigzag edges. The central section of the 
constriction operates as a single-level quantum dot, as the current flow towards the adjunct sections 
(rotated by 60°) is strongly suppressed due to mismatched valley polarization, although each section 
in isolation shows maximal quantum value of the conductance Go = 2e 2 /h. We further show, that 
the trapping mechanism is insensitive to the details of constriction geometry, except from the case 
when widths of the two neighboring sections are equal. The relation with earlier studies of electron 
transport through symmetric and asymmetric kinks with zigzag edges is also established. 

PACS numbers: 73.63.-b, 73.22.Pr, 81.07.Ta 



I. INTRODUCTION 

Soon after the breakthrough in its fabrication [TJ [2], 
an atomically-thin carbon monolayer (graphene) has at- 
tracted intense experimental and theoretical attention. 
An unusual band structure of graphene leads to exotic 
electronic properties [3], which makes possible either to 
create devices that have no analogue in silicon-based elec- 
tronics [4 , or to test various predictions of relativis- 
tic quantum mechanics in a condensed-matter system 
[5]. Additionally graphene's true two-dimensional nature 
combined with high carrier mobility makes it a promis- 
ing base material for studying low dimensional systems, 
such as quantum wires realized as graphene nanoribbons 
[SI H [1 ED] , or quantum dots [TTJ ECU EES] • Just to give 
a recent example, the exotic features of quantum chaotic 
behavior, characteristic for massless spin-1/2 fermions 
confined in a quantum dot [T3], have been found both 
in an experiment [T^] and computer simulations [T5] . 

Theoretical research on graphene nanostructures have 
started much earlier [TH EEI1 HH EEFJ but speed up after 
it was realized, that such systems are promising building 
blocks for a solid-state quantum computer [2QJI2T]. In at- 
tempt to operate on a solid-state qubit |22] in graphene, 
one need to deal with an obstacle that electrons occur in 
two degenerate families, corresponding to the presence of 
two different valleys in the band structure. Trauzettel et 
al. [20] propose to solve this problem by using the insulat- 
ing graphene nanoribbon with armchair edges, for which 
the valley degeneracy is lifted up by the boundary condi- 
tion [23] . Then, applying the gate voltage to a finite sec- 
tion of the ribbon, one forms the quantum dot which has 
a single relevant electronic level (with a spin-only degen- 
eracy) for a considerably wide interval of electron Fermi 
energies, and thus is called a single-level quantum dot 
(SQD). As a fabrication of perfect armchair nanoribbon 
seems to be difficult due to the edge instability [531 US] 
an alternative approaches, employing simultaneously ex- 
ternal magnetic field and the mass confinement to break 
valley degeneracy in graphene rings [S0J H7] , disks [25] , or 



antidots [121 [30], have been discussed. The operational 
simplicity of the device [2UJ , offering a fully-electrostatic 
control, has inspired designing its counterpart based on 
a nanoribbon with predominantly zigzag edges |31] , sim- 
ilar to these already obtained in well-controlled fabrica- 
tion processes [7j [8j [9] . An additional motivation to focus 
on ribbons with zigzag edges comes from the two theoret- 
ical findings: (i) local-density approximation results [32"] 
show that the electronic structure narrow graphene con- 
strictions can be well described by a simple tight binding 
model, and (ii) analytical discussion of the tight-binding 
equations for a semi-infinite honeycomb lattice [33 shows 
that the zigzag boundary condition applies gcncrically 
(at low energies) to arbitrary lattice termination, except 
from the case we have a perfect armchair edge. 

The device [3JJ, however, still needs sections of insu- 
lating ribbon with armchair edges, attached serially to 
both sides of the central section with zigzag edges (op- 
erating as SQD) to suppress the outward current flow. 
In this paper, we demonstrate numerically that a similar 
Z-shaped constriction, in which each of the three sec- 
tions has zigzag edges (see Fig. [ija) ) and carries a fully 
conducting channel, may operate as highly-effective SQD 
due to mismatched valley polarization between the neigh- 
boring sections (rotated by 60°). The operation of such 
a double-kink device is rationalized by referring to the 
conductance spectrum of a single-kink device (Fig. [l|b)) 
which suppress the current [33], and by the discussion of 
level quantization for a finite section of the ribbon with 
zigzag edges [35] , 

The paper is organized as follows: In Sec.[ll]we briefly 
present the numerical method of the conductance calcu- 
lation in a tight-binding model of graphene, and discuss 
its relation to the effective Dirac equation. Then, in Sec. 
|III[ we recall the idea of valley polarization in constric- 
tions with zigzag edges leading to the current suppression 
in the kink device, and calculate numerically the conduc- 
tance of kink devices with various geometries. Finally, in 
Sec. |IV| two kinks are attached serially and the resonance 
transmission via quantum-dot levels is discussed. 
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II. TIGHT BINDING MODEL FOR ELECTRON 
TRANSPORT IN GRAPHENE 

We start from the nearest-neighbor tight binding 
model taking into account the 2p z orbitals of carbon 
atoms [3], with the Hamiltonian 



tering matrix 



(i) 



The hopping matrix clement t^- = — t if the orbitals 
\i) and \j) are nearest neighbors on the honeycomb lat- 
tice (with t — 2.7 eV), otherwise r^- = 0. The elec- 
trostatic potential energy Vi = V(xi) varies only along 
the x-axis in the coordinate system of Fig. [I] It is cho- 
sen as Vj — Voo in the leads marked by shadow ar- 
eas (x < or x > L) or Vj = in the device area 
(0 < x < L). For a given Fermi energy E, the chemi- 
cal potential fij = E — V(xj) is equal to /loo = E — V^, 
or to Ho = E, respectively. Throughout the paper, we 
analyze the system conductance as a function of uo at 
Voo fixed such that fi^ corresponds to large number of 
propagating modes (the heavily-doped leads limit). 

In the limit of zero bias voltage, phase-coherent trans- 
port properties of a noninteracting system such as de- 
scribed by the Hamiltonian are encoded in the scat- 
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FIG. 1: The systems studied numerically in the paper, (a) 
Quantum-dot device made of zigzag nanoribbon of the width 
Wd and the length Ld, attached to similar ribbons of the 
width Wi rotated by 60°. (b) The kink device formed by two 
zigzag ribbons of widths Wl and Wr. Each of the systems 
is connected to heavily-doped graphene leads (shadow areas) 
displaced by the distance L along the a;-axis. 
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which contains the transmission t (t') and reflection r 
(r') amplitudes for charge carries incident from the left 
(right) lead, respectively. The conductance is determined 
by the Landauer-Buttiker formula 



G = G Trtt^ = 2 * r J2 T "> 



(3) 



where Go = 2e 2 /h is the conductance quantum, and T n 
is the transmission probability for the n-th normal mode. 
Apart from the simplest cases, when different transverse 
modes are not mixed by the transport |23[ I37j and the an- 
alytical solutions are available, one need to calculate the 
transmission matrix numerically. This can be achieved 
in two steps: First, we find the propagating modes in 



the leads ^ 



' Tn,- 



(i) 

where ?/>„ ± denotes the n-th 



incoming/outgoing mode in the i-th lead, with their self- 
energies £n,±- Then, the scattering matrix ([2]) is obtained 
from the Lee-Fisher relation 



S = -l-i^^V^- 



1 



H-E 



(4) 



The self-energy term £ = T"i'dia,g(^ ny+ )^''P' represents 
the leads, with V the coupling matrix of the leads to 
the contact region. The matrix v — diag(— 2Im£„ i+ ) con- 
tains normalization factors proportional to the propagat- 
ing velocities for the modes. The direct matrix inversion 
in Eq. Q is usually avoided by finding the S matrix 
via the recursive Green's function algorithm, also avail- 
able in a version for a multi-terminal geometry |38] . For 
the two-terminal geometries considered here, and for the 
number of lattice sites < 10 5 , the direct matrix inversion 
with standard numerical routines is also very effective. 
According to Eq. |4), the transmission matrix t depends 
on the properties of both constituents of the device under 
consideration: the leads and sample area. For the case of 
a weakly-doped graphene samples however, it was shown 
[3"5] that the transport of Dirac fermions is insensitive to 
the details of the leads, provided they carry a sufficiently 
large number of propagating modes. 

The method of the conductance calculation in a tight 
binding model of graphene, described briefly above, rep- 
resents an adaptation — to the honeycomb lattice — of the 
method developed by Ando for a square lattice jJUj. 
Albeit technically similar, from the fundamental point 
of view the two methods represents quite different ap- 
proaches to mesoscopic physics. Ando's approach starts 
from a discretization of the Schrodinger equation in order 
to solve the scattering problem for the geometries which 
are not tractable in a continuous limit. Here, we deal 
directly with the microscopic model of graphene, allow- 
ing one to extent the discussion on the situation when its 
effective model for low-energy excitations, given by the 
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FIG. 2: Schematic illustration of the mechanism of cur- 
rent suppression in the kink device. Top panel: Hexago- 
nal first Brillouin zone of graphene. The valleys centered at 
two inequivalent Dirac points, labeled K and K' , are shown 
schematically. Middle panel: Dispersion relation for the rib- 
bon with zigzag edges. For the lowest mode, there is one- 
to-one correspondence between the direction of propagation 
(indicated by arrows) and the valley isospin. The polarity 
switches (i) when changes the sign of chemical potential ft, 
or (ii) when rotates the constriction axis by 60°. Bottom 
panel: First Brillouin zone at the two sides of the kink de- 
vice. Solid (dashed) arrows indicate the direction of propaga- 
tion for K (K') valley. As corresponding points in the Bril- 
louin zone propagate in opposite directions at the two sides 
of the kink, electron cannot pass the device unless additional 
(microscopic) intervalley scattering processes are present. 

massless Dirac equation [3J, breaks down due to scatter- 
ing the carriers between valleys |41J . A finite-difference 
method for solving the Dirac equation on a square lattice 
was also recently developed )42] , 

III. THE KINK DEVICE WITH ZIGZAG EDGES 
A. Valley polarization and the electron transport 

To understand the physics ruling the conductance of 
devices shown in Fig. [T] one needs first to recall basic 
facts concerning electron transport through the constric- 



tion with parallel zigzag edges: a building block of each 
device considered. Such a constriction, also known as 
the valley filter |21j . was shown to produce, upon bal- 
listic injection of current, the noncquilibrium valley po- 
larization in a sheet of graphene attached. Motivated by 
the related analytical result for generic boundary condi- 
tions of graphene flakes {33j we have shown numerically 
[34,, that valley polarization is also produced by con- 
strictions with other edges, apart from the perfect arm- 
chair edges. These observations are rationalized as fol- 
lows: For a constriction long enough, the transport be- 
come one-dimensional, and crystallographic orientation 
of edges determines the direction of propagation in the 
first Brillouin zone (see Fig. |2j. As the lowest propagat- 
ing mode in such a nanoribbon lack the twofold degener- 
acy of higher modes (and may be arbitrary close to one 
of the Dirac points K or K ') the crystallographic orien- 
tation also determines the valley polarization of current 
passing the constriction (which is opposite for the con- 
ductance then for the valence band). 

Now, the two quite different devices consisting of two 
valley filters (of the opposite polarity) in series, can be 
constructed. By simple applying a step-like potential 
profile to the straight ribbon with zigzag edges, one ob- 
tains the electrostatically-controlled valley valve [Si] . Al- 
ternatively, one can build the kink device [34 43 . In both 
cases, the current suppression due to mismatched valley 
polarization produced by two constituents of the device, 
is expected. However, a bit more detailed discussion of 
the current-suppression mechanism on the example of the 
kink device (see Fig. [2j bottom panel) unveils its striking 
feature: Strictly speaking, a one-to-one correspondence 
between the direction of propagation and the valley in- 
dex, appearing in a nanoribbon, causes the electron can 
neither be transmitted nor reflected by the interface be- 
tween two filters of opposite polarity! The same obser- 
vation applies to the valley valve, for which the theoret- 
ical analysis of microscopic tight-binding equations |41j 
shows, that the intervalley scattering processes lead to si- 
nusoidal conductance oscillations around the mean value 
G = 0.5 Go when rotating the interface line with respect 
to the ribbon edge. The valve conductance also depends 
on its width, and is significantly different for the cases 
when zigzag and anti-zigzag ribbons are used, illustrat- 
ing the microscopic nature of an electron transport. For 
the kink device, an analytical solution is missing, and the 
existing numerical results are over vie wed briefly below. 



B. Conductance of the kink device 

The scattering problem for devices similar to that 
shown in Fig. [ljb) have been studied independently by 
several authors. The perfectly symmetric kink, formed 
by two semi-infinite nanoribbons with zigzag edges each 
of which is carrying a single propagating mode, shows 
irregular conductance oscillations with varying chemical 
potential [55] ■ The oscillations cover the full range of 
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FIG. 3: Conductance of the kink device as a function of the 
chemical potential at fixed L = 100\/3a (a) and as a function 
of the constriction length at fixed fio = A/2 (b). Inset in 
panel (a) shows a magnified region between the Dirac point 
and the first conductance step. In both panels, the width of 
left ribbon forming the kink is fixed at Wl = 40v3a, whereas 
the width of right ribbon is varied between the curves. 



splitting A = \\fZ-KTajWL = 7tt/80 and vl = 32 prop- 
agating modes for fx^ = Voo — t/2. The right ribbon 
width is varied as W R /(V3 a) = 10, 20, and 40. 

The operation of the kink device is demonstrated in 
Fig. [3] First, we took a sample area of the length 
L = lOQy/Sa and plot, in Fig.^a), the conductance as a 
function of the chemical potential. All three devices with 
different Wr show some conductance suppression at the 
first plateau \hq\ < A. The conductance of the symmet- 
ric kink (Wr = W L ) is still of the order of G ~ 0.1 G , 
but the asymmetric kinks shows much stronger suppres- 
sion at low energies (see the inset). For an additional 
illustration we present, in Fig. [3^b), the conductance as 
a function of L at fixed /io = A/2. For L = 0, we have 
G/Gq ~ min(j/£, vr) — vr, with the number of propagat- 
ing modes in the right arm vr = 32 or 8 (for Wr — Wl 
or Wl/4, respectively). For larger L, G first decrease 
exponentially, in order to saturate for L > \^3Wl/2 (the 
length above which the role of evanescent modes becomes 
negligible) in the symmetric case. For Wr = Wl/4, G 
continues to decay and approaches the value ~ 10~ 3 Go 
for longest examined systems, confirming the claim of 
Ref. |34j . that the asymmetric kink device in graphene 
blocks the current very effectively at low dopings. These 
findings also coincides with earlier results for Aharonov- 
Bohm quantum rings [26, 45 for which intervalley scat- 
tering rate was shown to be related to the presence of an 
asymmetry or edge irregularities. 

Away from the first plateau, each device shows con- 
ductance quantization that is essentially governed by the 
narrower arm width Wr. In the symmetric case, how- 
ever, we have the steps corresponding to even multiplici- 
ties of Go , whereas in the asymmetric case odd multiplic- 
ities appear. This can be summarize in an approximating 
formula for the kink conductance 



^ G ^ Go, and the upper limit is approached when the 
resonances with quasi-bond states, localized at the kink 
symmetry axis, occur. The system of tight-binding equa- 
tions describing the symmetric kink, however, was found 
to be numerically stable only for ribbons of a moderate 
width Wl — Wr < lO-s/3 a (in units of the lattice spacing 
a). In the case of an asymmetric kink with Wl/Wr = 2 
attached to heavily-doped graphene leads [31] we have 
G Go at the first conductance step, and the resonances 
do not appear. On the other hand, the kink-like system 
formed by joining the two ribbons rotated by 120° shows 
almost a perfect transmission, as the mismatched-valley 
polarization does not appear in such a case [44] . 

Here, we consider a slightly different geometry then 
studied in Refs. [331 S3]- Namely, the kink device is at- 
tached to heavily-doped graphene leads with armchair 
edges. This guarantees the reflection symmetry of the 
system when Wl = Wr, and eliminates the problem 
with numerical stability occurring for wider ribbons in 
the setup of Ref. [33]. We have fixed the left ribbon 
width at Wl = 40a/3o, corresponding to the subband 



G/G « mh\{v L ^R) - 8 VL , VR , (5) 

where the Kronecker delta S ULtUR accounts for mis- 
matched valley polarization of the lowest propagating 
modes in two kink arms, which affect the total conduc- 
tance also at higher dopings. 



IV. BOUND STATES OF A DOUBLE KINK 

The conductance spectra of Z-shaped nanoconstriction 
of Fig. [TJa) are presented in Fig. |4j We consider three 
different devices, all built by joining finite sections of two 
nanoribbons with zigzag edges: one of the 15v3a width, 
and the other of the 20-^/3 a width. In the first case, we 
arrange the building blocks in a double-kink setup, such 
that the central section is wider then the peripheral sec- 
tions (Wd — 20v3a, Wi = 15V3a). Next, we consider 
the setup of a uniform width (Wd — W\ — 15\/3 a). Fi- 
nally, we took the peripheral sections wider than the cen- 
tral one (W D = l5V3a, Wi = 20V3 a). For each setup, 
the total sample area length is fixed at L — 160 a. The 
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FIG. 4: Conductance of the double-kink device for different 
arrangements of its ribbon-like sections, (a) Narrow-wide- 
narrow, (b) uniformly- wide, and (c) wide-narrow-wide setup. 
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FIG. 5: Quantum-dot energy levels of the double-kink device. 
Lines: dispersion relation of a nanoribbon with zigzag edges. 
Datapoints in the inset: conductance maxima of a double kink 
with W D = 15V3a, Wi = 4W D /3, and L D = 110a. 



results show, either for the narrow- wide-narrow or for the 
wide-narrow- wide setup (see Figs. Qa) and[4jc), respec- 
tively), that G w at the first conductance plateau ex- 
cept from the narrow peaks, corresponding to resonances 
with the quantum-dot states of the central section. Only 
for a uniformly- wide device some finite intervals of 
for which G < G , appear in Fig.Qb). 

These findings are related to the operation of the kink 



device studied in previous section as follows: The asym- 
metric kink blocks the current at low dopings, so two 
identical devices of its kind in series function as an 
electrostatically-controlled quantum dot (with an elec- 
tron trapped between the kinks) regardless we join to- 
gether wider, or narrower of the kink arms. The two 
symmetric kinks in series may also trap an electron acci- 
dentally, when the eigenstate of the central section lies in 
an energy interval, for which the single-kink conductance 
is low enough. 

In Fig.[5j we plot the positions of the conductance max- 
ima extracted from Fig.Qc), together with the dispersion 
relation for an infinitely long nanoribbon of W = 15\/3& 
width. The fc-coordinate of each maximum is related to 
its number in sequence j as kj = — j)/Lr>, with 
j* determined via E(K) = A/2. (Notice that the first 
conductance maxima are blurred at finite machine pre- 
cision due to the dispersionless character of the lowest 
subband). Strictly speaking, each bound state of the 
central ribbon section of Ld length is given by a unique 
quantum superposition of propagating waves from K and 
K' valleys [35 , the dispersion relation for one valley is, 
however, a mirror reflection of that for the second valley. 
The datapoints shown in the inset of Fig.[5]follow the dis- 
persion relation for the lowest subband of a nanoribbon 
with zigzag edges. Hence, for each peak in the conduc- 
tance spectrum of a double-kink device, the correspond- 
ing quantum-dot energy level is found. 

V. CONCLUSIONS 

In conclusion, we have studied numerically the elec- 
tron transport through Z-shaped nanoconstriction with 
zigzag edges, operating as a quantum dot in graphene. 
The device conductance, analyzed as a function of the 
chemical potential, exhibits a series of narrow resonance 
peaks. Each of the resonances is linked up to the energy 
level of the central section, when separated from the other 
parts of the device. An electron-trapping mechanism is 
discussed by analyzing the transport through a basic de- 
vice building block: the kink with zigzag edges. A mod- 
erate current suppression is observed for the symmetric 
kink, whereas the asymmetric kink was found to reflect 
electrons almost perfectly. The role of a mismatched val- 
ley polarization on both sides of the kink is stressed. 

Similar bound states were recently found in the sim- 
ulation of transport through 5-shaped nanoribbon [15] 
with irregular edges. However, such a system contains 
finite sections of a nanoribbon with armchair edges, so 
the nature of the electron-trapping mechanism is not as 
clear as for the simpler system considered here. 
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